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Abstract. The alloy Ca2- a; Sr a ;Ru04 exhibits a complex phase diagram with peculiar magnetic metallic 
phases. In this paper some aspects of this alloy are discussed based on a mean field theory for an effective 
Kugel-Khomskii model of localized orbital and spin degrees of freedom. This model results from an orbital 
selective Mott transition which in the three-band system localized two orbitals while leaving the third 
one itinerant. Special attention is given to the region around a structure quantum phase transition at 
x ~ 0.5 where the crystal lattice changes from tetragonal to orthorhombic symmetry while leaving the 
system metallic. This transition yields, a change from ferromagnetic to antiferromagnetic spin correlations. 
The complete mean field phase diagram for this transition is given including orbital and spin order. The 
anisotropy of spin susceptibility, a consequence of spin-orbit coupling and orbital correlation, is a tell-tale 
sign of one of these phases. In the predominantly antiferromagnetic phase we describe a metamagnetic 
transition in a magnetic field and show that coupling of the itinerant band to the localized degrees of 
freedom yields an anomalous longitudinal magnetoresistance transition. Both phenomena are connected 
with the evolution of the ferromagnetic and antiferromagnetic domains in the external magnetic field and 
agree qualitatively with the experimental findings. 

PACS. 7 5.30.-m,75.50.-y,7f.25.Tn 



1 Introduction 

Layered ruthenate compounds have gained an important 
position as a class of transition metal oxides displaying di- 
verse physical properties. The interest in these materials 
has been initiated by the discovery of unconventional su- 
perconductivity in the single-layered Sr 2 Ru04 which rep- 
resents most likely a realization of spin-triplet Cooper 
pairing, an analog of the A- phase of superfiuid 3 He.0 
12] Since the mechanism for pairing is probably of mag- 
netic origin, the magnetic properties of related systems 
have been studied extensively. Increasing the number of 
RuC>2-layers, following through the Ruddelson-Popper se- 
ries Sr n+ iRu„03„+i with n as the number of layers per 
unit cell, the system turns into a ferromagnet for n > 3 
■EI Peculiar metamagnetic behavior was reported for the 
double-layer system (n = 2) and was interpreted as a novel 
type of quantum critical behavior. 0] On the other hand, 
replacing Sr by the isoelectronic Ca in the single-layer 
compound yields a Mott-insulator with antiferromagnetic 
long-range order of localized spin-1 degrees of freedom. 
These contrasting properties of the single- layer compound 
motivated Nakatsuji and Maeno (NM) to investigate the 
continuous series of alloys Ca2- a; Sr a ;Ru04.[5] The interpo- 
lating phase diagram is complex. NM identified three dis- 
tinct ranges of x characterized by different crystal struc- 



tures as well as electronic and, in particular, magnetic 
properties. The present study is concerned with a certain 
range in the NM phase diagram of these alloys (FigUJ. 
Therefore it is in order to briefly review here the basic 
experimental results for C^-xSr^RuO^. [5] 

Starting from the metallic and superconducting sto- 
ichiometric compound S^RuCU (x — 2) the introduc- 
tion of the smaller ion Ca for Sr leads to a contraction 
of the crystal volume. This is achieved by the rotation of 
the Ru06-octahedra, which occurs for x < 1.5. [Sj Never- 
theless, up to dopings of x = 0.5 in Ca2- 3; Sr 2 ;Ru04 the 
crystal structure remains tetragonal. This doping range, 
called region III by NM, is characterized by an increase 
of the uniform susceptibility which eventually becomes 
Curie-like with a Curie constant corresponding to nearly 
free spin 5 = 1/2 degrees of freedom as x = 0.5 is ap- 
proached. In the doping range 0.2 < x < 0.5, the region 
II following NM, the system suffers a further crystal de- 
formation through tilting of the RuOe-octahedra around 
a non-symmetric axis in the basal plane, leading to an or- 
thorhombic crystal symmetry. 6 This results in modified 
magnetic properties: the Curie-like behavior is truncated, 
i.e. the susceptibility decreases after a pronounced maxi- 
mum at low temperature. This behavior suggests antifer- 
romagnetic correlations with a characteristic temperature 
scale which gradually increases, from T ~ OK at a; = 0.5 to 
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about 10 K as x approaches 0.2. In addition, a pronounced 
inplanc anisotropy in the spin susceptibility appears, in- 
dicating spin-orbit coupling effects with a two-fold sym- 
metry axis. For both region II and III the system remains 
metallic. Finally, there is a first-order phase transition at 
x « 0.2 to a phase (region I, i.e. < x < 0.2 ) with flat- 
tened RuOg-octahedra at low temperatures. This phase is 
continuously connected with the pure Ca2RuOg, which is 
a Mott-insulator and shows antiferromagnetic order with 
Tm ~ 100 - 150 K. The unusual sequence of changing 
properties upon Ca-Sr alloying has motivated a number 

Before going into details of different attempts to ex- 
plain the physics underlying these properties we discuss 
the basic electronic structure of the metallic Sr2Ru04, 
which is well understood experimentally as well as theore- 
tically. 7 8 : 9 This will be the basis for all following discus- 
sions of the alloy. The electronic band structure receives 
its basic character from the 4c£-i2 9 orbitals of the Ru-ion: 
{d yz ,d zx ,d xy } which disperse dominantly in each RuC>2- 
plane via the 7r-hybridization with the 0-2p-orbitals. It 
is easy to see that the d yz - and d^-orbitals have a very 
anisotropic dispersion yielding two quasi-one-dimensional 
bands, that hybridize to form one electron-like (/3-band) 
and one hole-like (ev-band) Fermi surface. The quasi-one- 
dimcnsional nature of these bands produce strong nesting 
effects which gives rise to enhanced incommensurate spin 
fluctuations around a wave vector Q w (27r/3a, 27r/3a, 0) 
|10U11II12| as observed in neutron scattering (a being the 
inplane lattice constant) . |13| The 7r-hybridization with the 
intermediate O 2p-orbitals in both inplane directions gen- 
erates for the d^-orbital a wider genuinely two-dimensional 
7-band which has van Hove singularities rather close to 
the (electron-like) Fermi surface. In this way the band 
structure known from de Haas- van Alphen [7] and ARPES 
measurements [E] is qualitatively well understood. 

The doping of Ca for Sr introduces a gradual rota- 
tion of the RuOg-octahedra around the z-axis and yields 
a modification in the band structure, which the analysis 
of the emerging magnetic properties in region III is based 
on.[S] Both groups, Nomura and Yamada (NY) [TO], and 
Fang and Terakura (FT)J3], emphasize that octahedra 
rotation tends to narrow particularly the 7-band affecting 
the band structure in a way that the 7-band van Hove 
singularity becomes a dominant feature for the electronic 
properties. This leads to an increase of the uniform spin 
susceptibility in the region III via the increase in den- 
sity of states and the Stoner enhancement. On the other 
hand, octahedra tilting, as occurring in region II, tends 
to narrow the a-/3-bands and would shift the nesting vec- 
tor Q = (7r, 7T, 0) pushing the system more towards an 
antiferromagnetic instability. While this scenario leads to 
a system with metallic behavior and an enhanced uni- 
form susceptibility in region III, it is difficult to gener- 
ate the Curie-behavior displaying apparently one localized 
S = 1/2 degrees of freedom per Ru-ion. This spin size is 
even more surprising in view of the fact that the fully 
localized 4t2 9 -orbitals with four electrons on the Ru A+ - 
ion would by Hund's rule favor a spin-1 configuration, as 
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Fig. 1. Basic schematic phase diagram of Ca2-a;Sr :l: Ru04 at 
low temperatures (following Nakatsuji and Maeno|5]). The re- 
gion III is metallic, paramagnetic and has tetragonal crystal 
symmetry. The region II is metallic, at low temperature with 
presumably antiferromagnetic correlation and has orthorhom- 
bic symmetry. In both regions the susceptibility over a certain 
range of temperature is Curie-like with a spin S = 1/2. The 
region I is a Mott insulator with antiferromagnetic long range 
order of a spin 5 = 1. 



observed in Ca2Ru04. Of course, complete localization is 
ruled out for the system being metallic. 

In contrast to the NY and FT models based on three 
modified itinerant electronic bands, Anisimov et al. fo- 
cused their attention on the emerging spin S = 1/2 in 
the susceptibility around x = 0.5. |15U16| In their analysis 
the important feature of the RuOg-octahedra rotation is 
a progressive narrowing of the bands with increased Ca- 
doping. Due to the different character of the three bands 
(the parity with respect to reflection at the basal plane is 
opposite for the a/3-orbitals and 7-orbital), electron cor- 
relation will then drive a Mott transition which is orbital 
selective. Since the a-/3-bands, derived from quasi-one- 
dimcnsional bands, have in Sr2Ru04 approximately half 
the width of the genuinely two-dimensional 7-band, cor- 
relation would be more effective to localize their carriers. 
Using LDA augmented by dynamical mean field theory 
(DMFT) and non-crossing approximation (NCA) Anisi- 
mov et al. demonstrated that the band narrowing can 
indeed lead to a redistribution of the four electrons per 
Ru-ion and Mott transition for the a- and /3-band.[T51 
I16| The two orbitals d yz and d zx absorb 3 electrons while 
the 7-band with one electron, although half-filled, remains 
metallic. The 3 electrons on the two orbitals d yz and d zx 
generate a localized spin 1/2 and an orbital isospin 1/2 
degrees of freedom. This change of electronic properties is 
assumed to occur upon doping and is essentially reached 
when x — 0.5 is approached. For the range 0.2 < x < 0.5 
the magnetic properties are dominated by the combina- 
tion of orbital and spin degree of freedom. While this pic- 
ture results in a localized spin 1/2 per site, as observed 
in the spin susceptibility and leads to a slight elonga- 
tion of RuOg-octahedra (~ 1 — 2%) as the orbital oc- 
cupation is rearranged fB|, it seems in contradiction with 
other measurements. Recent optical spectroscopy data do 
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not observe as drastic a change of carrier concentration 
in the range 0.5 < x < 0.9 as one would anticipate from 
an orbital-selective Mott transition.^] Furthermore, an 
apparent inconsistency in the orbital structure has been 
reported based on neutron scattering data.JU| Recently, 
Licbsch argued based on a modified DMFT approach that 
orbital-selective Mott transition is unlikely to occur for the 
Ad — £2g-orbitals.[2D] However, Koga et al. have very re- 
cently demonstrated that under rather general conditions 
an orbital-selective Mott transition in a two-band Hub- 
bard model is possible, contradicting Liebsch's argument 

ED 

Despite the mentioned reservation towards the orbital- 
selective Mott transition we will use this scenario in this 
paper as a working hypothesis and investigate some possi- 
ble consequences mainly for the region II of the NM phase 
diagram. We find properties of the model which are in 
surprisingly good qualitative agreement with important 
features experimentally observed. Assuming three local- 
ized electrons in the d y2 -d zx -orbitals we discuss the ef- 
fective Kugel-Khomskii-type Hamiltonian describing the 
corresponding localized degrees of freedom, spin and or- 
bital isospin, within a mean field theory. Including doping 
disorder into the model gives an important clue for the 
comparison with experimental findings. In particular, it al- 
lows to discuss the observed metamagnetic transition and 
the unusual magnetoresistance behavior in region II. We 
will conclude that the picture of localized d^-d^-orbitals 
gives a good description of the magnetic properties of the 
region II which is probably difficult to explain with only 
itinerant electrons in all three orbitals. 



Hilbert space represented by an isospin. We define, there- 
fore, isospin operators / which act in the following way 

i z \±) = ±\\±), / + i-> = n+) = n- 

(i) 

Hence, on every site the localized state consists of a set of 
four product states 



2 Effective Kugel-Khomskii Model 

For doping concentrations x in region II and at the bound- 
ary of III the G? yz -G? z:r -orbitals are considered to be local- 
ized, according to Anisimov et fflL |15lllf>| We now intro- 
duce an effective model for these localized orbitals, and 
will neglect the still itinerant 7-band (which will play a 
role again in a later stage). This is obviously a drastic sim- 
plification in view of the influence the 7-band may have 
on the two localized orbitals. There are two basic types of 
coupling: (1) RKKY-interaction and (2) double exchange. 
The former tends towards an antiferromagnetic spin cor- 
relation, while the latter, in contrast, would prefer ferro- 
magnetism. Despite the fact that the 7-band is half-filled, 
the band structure is far from particle-hole symmetric, so 
that no dangerous singular features in the RKKY coupling 
would appear. Thus it is unclear which trend would win. 
We shelve this problem for the following discussion, ac- 
knowledging the short-coming of in our simplified model. 
The considerably more complex complete model will be a 
matter of future study. 

The three localized electrons in the two orbitals (d zy , d zx ) 
create four different configurations, consisting of a spin 
1/2, I |) and | j), and an orbital degree of freedom, |+) 
and |— ) describing the singly occupied d zx and d yz , respec- 
tively. Analogous to the spin also the orbital degree of free- 
dom corresponds to a two-dimensional SU(2)-symmetric 



{ITU IT-), U+>, II-)}. 



(2) 



The fundamental microscopic model providing the dynam- 
ics for these degrees of freedom is an extended two-orbital 
Hubbard model (neglecting the 7-band) where we restrict 
ourselves to nearest-neighbor hopping and onsite interac- 
tion for the intra- and inter-orbital Coulomb repulsion, U 
and U 1 , respectively, and the Hund's rule coupling Jh- 



T~La,0 — t { c i+a y ,yz,s c i,yz,s + c i+a m ,yx,s C h^>s + h.C.} 
— A* y ] c i,v,s C i,v,s 

,zx 

i v % 

+ Jh C \,yz,s C \,zx,s' C i,zx,s C i,yz,s' ' 

i,s,s' 

(3) 

where ; (c m .i iS ) creates (annihilates) an electron on 
site i with orbital index v (= yz, zx) and spin s (n.^.s = 

c i,u,8 C h",s> n h" = n wT + n i,v,l> a = ( a x,a y ) = (1,0) or 
(0, 1) basis lattice vector). The effective nearest-neighbor 
Hamiltonian in terms of spin and isospin has the form of a 
Kugel-Khomskii model and can be derived within second 
order perturbation in t/U, 



/.a 

with the coefficients 

3a 2 + 1 



.4 
B 
C 
D 



(3a - l)(a+ 1) 

-(1-a) 2 
(3a 2 + l)(3a - l)(a + 1) 
5 — 3a 



4(3a- 1) 
1 



(5 - 3a) (3a - 1) 



(3a-l)(a + l)^ 2 
2(3a 2 + 1) 
3a — 1 



2(5 - 3a 



\a 2 x - a 2 ) 



(4) 

(5) 
(6) 
(7) 
(8) 

(9) 
(10) 



4 



Manfred Sigrist, Matthias Troyer: Orbital and spin correlations in Ca2- :E Sr a ;Ru04: A mean field study 



where ry a and rj' a have opposite sign for the x- and y-axis 
bonds. We impose the (approximatively valid) relation 
U = U' + 2J H with a = U'/U and J = 4t 2 /U. In or- 
der to obtain a valid approximation a should lie between 
1/3 and 1, and not too close to these boundary values. We 
choose throughout this paper a = 3/4 as a representative 
value whenever we do a concrete analysis of this model. 
We emphasize that while the spin has complete SU(2) 
symmetry, the isospin has only Ising-like interactions, i.e. 
quantum fluctuations are suppressed for the latter. 

The Hamiltonian (TJJ describes the system with tetrag- 
onal crystal symmetry corresponding to regime III in the 
experimental phase diagram close to x = 0.5. Now we add 
terms which describe the orthorhombic distortion due to 
the tilting of the octahedra with the rotation axis in the 
basal plane. We do not speculate about the origin of this 
lattice deformation, which we consider to lie outside our 
model. Following Ref. ^1 we introduce these with the 
strains e\ — e xx — e yy corresponding to a rotation axis like 
[100] and £2 = e xy corresponding to the axis [110]. These 
distortions lift the orbital degeneracy of (d yz ,d zx ) and, 
hence, can be viewed as an effective "uniform field" polar- 
izing the isospins in a certain direction. From symmetry 
considerations it follows that 



H dist = Y^K^t + ] 



(11) 



where are coupling constants (see the Appendix). 

The first term is a longitudinal field while the second term 
corresponds to a transverse field, introducing quantum 
fluctuations. In connection with the Ising Hamiltonian 10} 
this second term provides the possibility of a continuous 
quantum phase transition. [22] The lattice distortion in the 
regime II has a tilt axis lying between [100] and [110], so 
that both strains are turned on, if the border at x = 0.5 
between the two phases is crossed. 



3 Mean field approximation 

We now discuss the basic properties of this model and 
its phase diagram within a mean field approximation for 
both the spin and isospin degree of freedom. The Hamilto- 
nian on the square lattice possesses a bipartite structure, 
which is the basis of our mean field decoupling. We intro- 
duce different mean fields for the corresponding A- and 
£>-sublattice. 



(if) 



t A i eA 

t B i G B 



and (Sf) 



sa i G A 
s B i G B 



(12) 



The mean field calculation is straightforward, whereby the 
self-consistent equations are solved numerically. We find 
the following properties. 



3.1 Tetragonal system 

First we consider the situation in the absence of any or- 
thorhombic strains, e\ = £2 = 0. For the choice of param- 
eters (a = 3/4) we find that the coupling of the isospins, 



represented by JC, has a higher energy scale than the 
one for the spins, leading to a higher mean field transi- 
tion temperature. The parameter C in Ea. l|4l7|l is positive, 
favoring a staggered (Ising-antiferro-) orbital order. This 
orbital order implies ferromagnetic exchange between the 
spins. Note that this gives rise to enhanced uniform spin 
susceptibility as observed in experiment. The transition 
temperature to ferromagnetic order is low compared to 
the energy scale of the orbital order, so that it is natural to 
expect over a wide range of temperature a Curie-like sus- 
ceptibility originating from almost independent spin 1/2 
degrees of freedom, consistent with the experimental ob- 
servation. 



3.2 Orthorhombic distortion 

Now the uniaxial strains £1 and £2 are turned on. Since in 
regime II the crystal symmetry is determined by the rota- 
tion of the octahedra around an inplane axis which does 
not coincide with a crystal symmetry axis, both strains, 
£1 and £2 appear. For our calculations we assume that 
Kei = Q.7bKt2 — Ke which we use as a control parame- 
ter for our phase diagram. Let us first consider the orbital 
part only, ignoring the spin part. We observe a competi- 
tion between the staggered order and the uniform align- 
ment induced by the strain field. The transverse isospin 
component I x induced by the strain £2 reduces the transi- 
tion temperature, driving the system towards a quantum 
phase transition. This is included in our mean field ap- 
proach with the following self-consistence equations. 

2JCt B + K 1 e 1 



t A = 

2y<{2JCt B + Kjei)^ + K%e$) 

x tanh {b^I {2JCt B + K^f + K$f% 

The isospin component t x = (I x ) is obtained as 



(13) 



txA — 



2 v / {2JCt B +K 1 e 1 ) 2 +K$4 



(14) 



x tanh l3J(2JCt B + K^) 2 + K^i) 



for the A-sublattice site. The analogous expressions for 
the B-sublattice site are obtained by exchanging t B and 
tA in these equations. The ^-component is only finite, 
if the orthorhombic strain £ 2 is present. The staggered 
orbital moment t s = {tA — t B )/2. Increasing the strains 
leads to a gradual diminishing of t s and the corresponding 
transition temperature until it would vanish at a quantum 
critical point. In turn the mean fields (t X A + t xB )/2 and 
^0 = {tA + £b)/2 increase continuously favoring AF spin 
correlation resulting in the phase diagram shown in Fig. [2] 
The decrease of t s leads at the same time to a weakening of 
the ferromagnetic correlation and the gradual increase of a 
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Fig. 2. Mean field phase diagram as a function of tempera- 
ture T and strain Ke. Thin lines denote second order phase 
transitions and the thick line a first order phase transition. 

uniform isospin correlation gives rise to a competing anti- 
ferromagnetic spin correlation. Indeed the continuous evo- 
lution of t s is truncated by the onset of antiferromagnetic 
spin order, before the anticipated quantum suppression of 
the staggered isospin order is completed. This spin corre- 
lation is favored by the strain driven ferro-orbital correla- 
tion and has a higher energy scale (mean field transition 
temperature) than the ferromagnetic order. Consequently, 
the staggered isospin phase and the antiferromagnetic or- 
der are separated by a first order transition. There is even 
a reentrant behavior associated with this first order in- 
stability such that around a critical range of Ke » 0.5 J 
upon lowering temperature first a staggered isospin phase 
with enhanced tendency towards ferromagnetism would 
be reached and with a first order transition (thick line in 
Fig. [21 at lower temperature the antiferromagnetic phase 
would appear. Otherwise all phase transitions in the phase 
diagram of Fig. are second order (thin lines). 

The evolution of the system, under orthorhombic dis- 
tortion, from ferromagnetic towards antiferromagnetic be- 
havior displays the correct qualitative trend in comparison 
with the experimental situation. The fact that the temper- 
ature scale of the antiferromagnetic state is higher than 
that of the ferromagnetic is also in qualitative agreement 
with experimental data. The former is easily visible in ex- 
periment at a temperature T ~ 5 — 10K around x — 0.2, 
while the latter can brought into connection with a irre- 
versibility transition observed in the vicinity of x = 0.5 
below ~ 1^.(23 The irreversibility and the sensitivity to 
slow dynamics led to the interpretation as a cluster glass. 
Similar features could arise in an inhomogeneous ferro- 
magnetic phase. |23 

3.3 Spin-orbit coupling and anisotropic spin 
susceptibility 

For the Ru-ion spin-orbit coupling is not negligible. The 
microscopic formulation of spin-orbit coupling involves the 



whole set of i2g-orbitals including the itinerant 7-orbital. 30 
We remain within our present reduced model and restrict 
ourselves to a phenomenological approach where the effect 
of spin-orbit coupling enters via the g-tensor modifying 
the Zeeman term in the Hamiltonian. This term can be 
derived from symmetry considerations analogous to the 
coupling of the strain (Appendix A) . 

Hz = -VB ]>> H ' Si + 9i{If)(H x Sf - H y Sf) 

(15) 

+g 2 {If){H x Sf + H V S?)] 

where only the two isospin directions I z and I x have been 
included, ignoring I v which is not induced by strain nor 
interaction (fin - Bohr magneton). Two coupling constants 
gi/g and g 2 /g enter as phenomenological parameters. (I a ) 
denotes the average of the orbital configurations on all 
sites. Here we want to concentrate on the anisotropy of 
the susceptibility in the basal plane (x-y) which shows 
the most distinguished effect. 

We now calculate the spin-susceptibility on the back- 
ground of the staggered orbital order and the driven ferro- 
orbital correlation: 

X (T, 0) = X0 (T){1 + K (P) cos 20 + * (I x ) sin 20} (16) 
9 9 

where the angle corresponds to the orientation of the field 
in the basal plane relative to the x-axis. The isotropic 
susceptibility Xo(T) includes the spin correlation and ap- 
proaches the behavior of free spins for temperatures much 
higher than the Need temperature in our model. 

In the tetragonal phase where a staggered Ising orbital 
order is realized the average of the isospin vanishes and, 
consequently, there is no basal plane anisotropy in the 
spin susceptibility. This is fully consistent with the exper- 
iment at x = 0.5. On the other hand, the orthorhombic 
distortion which yields polarized orbitals in region II gives 
rise to anisotropy. Let us for simplicity now assume that 
the orthorhombic distortion is of the type e xy ^ and 
e IX — e yy = 0. Hence we find that only (I x ) 7^ 0. If the 
temperature is much higher than the characteristic energy 
scale for the orbital correlation, the isospin can be approx- 
imately considered as an independent degree of freedom in 
a driving field, i.e. (I x ) w K 2 e xy j '4fcsT '. Here we assume 
that the strain e xy only weakly depends on temperature 
and keep it constant. This is only valid far away from the 
phase boundary of the structure phase transition. Thus, 
the anisotropy of the spin susceptibility can be described 
by the ratio of the maximal to minimal value 

= g + g 2 (I x ) K 4gk B T + K 2 e xy g 2 
X(-tt/4) g-g 2 (I x ) ~ igk B T - K 2 e xy g 2 ( ' 

which tends to 1 for large temperatures and grows for de- 
creasing temperatures. Obviously the anisotropy increases 
too, if e xy becomes larger, which implies a growing anisotropy 
for samples with x — > 0, 2 in region II. For a qualita- 
tive view of the behavior we plot for different parameters 
K 2 e 2 (FigEJl- This behavior is close to the observed one 
by NM.E3 
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Fig. 3. Anisotropy of susceptibility: Ratio of the maximal and 
minimal susceptibility as function of temperature for gilg — 
0.1 and different levels of distortion K^2- 

4 Effects of an external magnetic field and 
alloy inhomogeneity 

We focus now on the phase in which the lattice distor- 
tion induces the FO-order and AF-spin order. Close to 
the first-order phase transition line the phase with a stag- 
gered orbital order and the ferromagnetic spin correlation 
forms a metastable states. This can be stabilized by an 
external magnetic field, since the antiferromagnetic phase 
has a lower uniform spin susceptibility. A discontinuous 
"metamagnetic" transition appears, since the magnetic 
field allows us to traverse the first-order phase transition 
line and to reach the phase of higher magnetization. This 
type of transition had been investigated in another con- 
text by Khomskii and Kugel [21] ■ In our mean field cal- 
culation we apply the magnetic field along the z-axis and 
introduce two additional mean fields s x a and s x b as or- 
der parameters of the antiferromagnetic order perpendic- 
ular to the field, in addition to the uniform component 
along the field. In Fig. 0] we show the mean field results 
of the magnetization curves for three temperatures and a 
scries of distortion strengths. For low temperature indeed 
a discontinuous transition between a low and high mag- 
netization state is found. Obviously, this transition would 
be hysteretic since it is first order. Obviously, the field- 
induced first-order transition leads to a rearrangement of 
the orbitals and should result in a structural deformation 
as well 



4.1 Metamagnetic transition 

The real material does not exhibit a discontinuous meta- 
magnetic transition even at low temperatures. The short- 
coming in our model lies in the assumption of homogene- 
ity. Obviously, the random alloy would rather have spa- 
tial fluctuations, for example, in the crystal deformation 
field Ke. A spatial variation of Ke would give rise to the 
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Fig. 4. Magnetization curves at three temperatures for differ- 
ent values of Ke/J. 
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appearance of both the antiferromagnetic and the stag- 
gered orbital phase with ferromagnetic spin fluctuations 
in form of domains. They are separated by rather sharp 
domain boundaries, if the randomness is sufficiently broad 
and smooth. |2S] We distinguish A- and i^-domains which 
denote the domains with antiferromagnetic and ferromag- 
netic spin correlation, respectively. We consider a system 
where the antiferromagnetic phase dominates and the F- 
domains form at most small islands, as would be likely 
the case for x « 0.2 based on our model. These .F-islands 
develop a magnetic moment at low temperature once the 
ferromagnetic correlation length reaches their spatial ex- 
tension. In turn the moments of different islands are corre- 
lated via the antiferromagnetically correlated background. 
This effective interaction may be frustrating so that un- 
der certain conditions at low enough temperature even 
glass-like behavior can occur. We will here, however, con- 
centrate on the behavior under a magnetic field. 

The external magnetic field H ex polarizes the mag- 
netic moments of both domains. Because the F-domains 
possess the larger susceptibility, their free energy density 
drops more rapidly with increasing field than that of the 
A-domains. Thus the F-domains grow at the expense of 
the A-domains. The initial linear response of the magneti- 
zation to the field turns nonlinear, if the domains change 
their size. In this way the metamagnetic transition, which 
was discontinuous and hysteretic in the homogeneous sys- 
tem, would be smeared, continuous and reversible for the 
inhomogeneous case.j2El The inflection point of the mag- 
netization as a function of H ex , would roughly correspond 
to the highest density of the domain boundaries through 
the system, because this yields the most rapid change of 
the domain sizes as a function of H ex . The metamagnetic 
transition lies, consequently, also rather close to the point 
where the domain boundaries percolate throughout the 
sample. 



4.2 Extended mean field approach 

The continuous metamagnetic behavior of the inhomoge- 
neous system can be very easily simulated within our mean 
field calculation. We modify Ea. HllI) to a random field cou- 
pling 



H 



dist 



E 



[K 



K 2 e 2i I? 



(18) 



so that for every site, in principle, the crystal field may be 
different. Since one has to consider certain realizations of 
randomness we can only deal with finite lattices. We use 
a lattice size of 100 x 100 to simulate the evolution of the 
domain distribution in a magnetic field. For this purpose 
the model was set up in the following way. In order to 
obtain a smoothly varying distribution we start with an 
uncorrelated uniform distribution of random numbers in 
the interval [—0.9, 2.7] Averaging the values five times with 
the values at the four nearest neighbor sites we finally 
arrive at a smooth distribution shown in Fig. El with mean 
value K e/ J = 0.9 and with a standard deviation AKe/j = 
0.2, correlated over several lattice sites. 




Fig. 5. Distribution of Ke/C used in the mean field calcula- 
tion. The mean value of Ke/C is 0.9 and the standard deviation 
0.2. Black corresponds to the minimum value of 0.2 and white 
to the maximum value of 1.6 



For the mean field calculation we start from an initial 
configuration that has staggered orbital order for the t z 
mean fields, uniform magnetic order for the s z mean fields 
in the direction of the magnetic field and staggered anti- 
ferromagnetic order for the s x mean fields perpendicular 
to the magnetic field. In order to find a zero-temperature 
mean field solution we then perform 10 6 updates to this 
mean field, solving the mean field equations each time for 
one randomly chosen spin, using a quantum operator li- 
brary developed by one of the authors • 

Since most of the K e values are inside the antiferro- 
magnetic phase we find only small islands of staggered or- 
bital order and ferromagnetism (Fig. EK). Increasing the 
fields we can observe that these ferromagnetic domains 
grow, and percolate at a "transition" to a ferromagnetic 
phase (Fig. Eb,c). This transition is accompanied by a 
smooth metamagnetic like magnetization curve, shown in 
Fig. [7^. The volume of the F-domains grow monotonically 
as shown in FigdD. 

We would like to mention here that the presence of a 
metamagnetic transition has been observed in region II 
28,22). Also a structural response has probably be seen 
in experiment |32j . 



4.3 Phenomenological model 

We now illustrate the result of our mean field simulation 
by means of a simple phenomenological model describing 
the distribution of the F- and A-domains. We define the 
free energy at a given temperature as 



W 



F(l;H ex ,T) = ^-{l-h? -V F (l)f F (H ex ,T) 
V A (l)f A (H ex ,T), 



(19) 
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Fig. 6. Distribution of ferromagnetic (black) and antiferro- 
magnetic (white) domains at a) a weak field h/J — 0.2 and 
b) close to the percolation threshold with h/J — 0.3 and c) at 
stronger fields h/J = 0.4. 
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Fig. 7. Zero temperature magnetization curve and volume 
fraction of the ferromagnetic domain of the inhomogeneous 
system with the distribution of Fig. 



where we denote the integrated area of the F- and A- 
domains by Vp(l) and Va(1) (note Vp + Va — V the total 
area), and I is a variable determining the domain distri- 
bution. It is chosen so that the energy expense by devi- 
ating from the zero-field equilibrium state (characterized 
by I = I q) follows the quadratic behavior with an "elastic 
constant" W. Increasing I > Iq leads to the enlarging of 
Vf(1) which is a monotonically growing function of I. The 
free energy densities in the two domains are denoted as 
/f and Ja and are a function of the field. 

The equilibrium value of I for a given external field 
H ex is obtained by minimizing the free energy. This leads 
to the equation 



W(l - h) + [f F (H ex ,T) - f A {H ex ,T)]v{l) = (20) 
with 
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dV F 



(I), 



(21) 



a measure for the density of domain boundaries. It is ob- 
vious that v(l) would be non-monotonic with a maximum 
at an intermediate value of I. The integral J dlv{l) run- 
ning from I = — oo to +00 gives the total volume V . This 
function characterizes the inhomogeneous system within 
this theory. Since fp < Ja in a finite field the solution of 
Ea. H4.3JI yields l(H ex ) as a monotonically growing function 
of H ex . The magnetization of the system is given by 



M 



dF 

dHprr 



-V F (l) 



df F 
dH„ 



(v - VHO) 



dH P 



(22) * 



= V F (l)M F + (V- V F {l))M A 



We use the following simple approximative forms for the 
free energy densities in a magnetic field: 




fa {Hex , T) — fo a 



M 2 
IU 0a 

XOa 



lncosh(xoaH ex /M 0a ) (23) 



with a = F,A and f 0a = / o (0) (/of = /oa)- Here M 0a 
is the saturation moment and Xao is the uniform spin 
susceptibility in the linear response regime at the given 
temperature T. Thus, the magnetization is given by 



M„ 



dfa 

dH e .T 



M 0a t&nh(xoaH ex /M 0a ) (24) 



Note that in this discussion we do not rely on the presence 
of long range order, neither FM nor AFM. The character- 
istic function v(l) depends on the specific realization of 
the alloy. We may take as a simple example 



/(!) = 



V 



L 



2arctan(Z) I 2 + L 2 



(25) 



for |/| < I and v{l) = otherwise and — I < Iq < 0. 
Ea. l|4.3|l can solved numerically. For a given set of pa- 
rameters we display the results in Fig. |SJ In Fig. the 
volume of the F-domains as a function of H ex is shown 
to grow monotonously and in Fig. [Hp) the magnetization 
exhibits a clear metamagnetic transition at the point of 
fastest change of the volume. 



4.4 Longitudinal Magnetoresistance 

We now would like to consider the effect of domains on 
the transport in the metallic 7-band which we had ne- 
glected so far. Naturally the conductance of the 7-band 
depends also on the correlation of the localized degrees of 
freedom of the a-/3-orbitals, since they interact via onsite 
couplings, such as Hund's rule coupling. First of all the 
two uniform phases A and F would have slightly different 
electrical conductance, a a and a F . For the inhomogeneous 
phase the domain boundaries also play a role. They are 



Fig. 8. Results of the phenomenological model: a) Volume of 
the F-domain; b) magnetization ; c) longitudinal magnetore- 
sistance. Parameters: W = 0.1 V, L = 2, lo = 2, I — 2.5, 
M 0F = M A = 1, xof = 0.3 and X oa = 0.02. 



restricted regions where orbital and spin degrees of free- 
dom have a rapid spatial dependence connecting the two 
kinds of domains, so that the scattering of electrons of 
the 7-band in the domain boundary region is enhanced. 
This readily understood from the point of view that the 
A- and F-domains in a magnetic field have different mag- 
netizations. Thus the rapid change of magnetization in the 
domain boundary constitutes a scattering potential. 

We would like to consider the problem of the longitu- 
dinal magnetoresistance, i.e. the resistance for a current 
which flows parallel to the applied external (inplanc) field. 
The change of the domain distribution is the key ingredi- 
ent determining the field dependence of the resistance. We 
address this problem again from the point of view of our 
simple phenomenological model. The overall conductance 
of the sample neglecting the domain boundaries may be 
roughly approximated by the averaged conductance 



V F a F + (V- V f )(ta 
V ' 



(26) 



if the conductances a F and a F have comparable magni- 
tude. Next we introduce the effect of the domain bound- 
aries. As additional scattering regions their influence de- 
pends on their density. Hence, the resistance contains two 
parts coming from the domains and the domain bound- 
aries: 

p. (27) 



1 

a 



Here p denotes the domain boundary part which we ap- 
proximate as proportional to the density of boundaries, 
corresponding to v. p = p^v. Now we can discuss the 
field dependence of p by inserting the solution of Eq. (|4.r?p , 
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l(H ex ). This leads to a sensible qualitative behavior of the 
longitudinal magnetoresistance as can be seen in FigJHt) 
showing Ap(H ex ) — p(H ex ) — p(0). There is a pronounced 
maximum at essentially the same position as the metam- 
agnetic transition occurs, both connected with the domain 
boundary percolation condition. It is also worth noting 
that the low-field limit leads to 

Ap(H ex )<xl-l cxf F (H ex ) - f A (H ex ) cx H 2 ex . (28) 

This low-field behavior and the overall field dependence 
compares well with the experimental data by NM. |27 , 28 
In the experiment a pronounced maximum is found at ap- 
proximately the field corresponding to the metamagnetic 
transition. This strong feature in the magnetoresistance 
disappears for temperature higher than the characteris- 
tic temperature of the antiferromagnetic spin correlation 
(~ 1QK). Naturally the domain formation is absent when 
there is no spin correlation (AFM order in our mean field 
treatment). 



4.5 Resistor network approach 

The issue of transport in such a system is naturally com- 
plex since it involves percolation of domains and domain 
boundaries. This aspect is not taken into account in our 
model where we the domain boundaries only considered as 
additional uncorrelated scatterers entering the resistance 
only through their density. The domain boundary, how- 
ever, should be considered as a barrier the electron would 
have to traverse. Hence, it would redirect the current flow 
which is difficult to include in our simple model. In or- 
der to demonstrate, however, that the basic properties of 
the electrical resistance are not changed by including this 
aspect, we consider a two-dimensional resistor network 
model where we can implement the progressive change of 
the conductivity and the presence of domain boundaries 
as barriers in a simple wav.|31| 

The network is a system of N x L knots on square 
lattice which coupled to their four nearest neighbors via 
resistors. The knots lie either in the A- or F-domain. Two 
neighboring knots in the same domain are connected via 
a resistor with the corresponding domain resistance Ra 
or Rp , respectively, while the resistor between neighbor- 
ing knots of different domains has a resistance Rdb which 
is considerably larger, corresponding to the barrier effect 
of the domain boundary. We assume a potential difference 
along L-direction of the network and apply Kirchhof 's law 
which yields a linear current voltage relation, i.e. the re- 
sistance. The domain distribution is generated by means 
of a real function Q(n x ,n y ) where n x and n y are integers 
denoting the knots. Introducing 

Q(n x ,n y ;l) = Q(n x ,n y ) + I (29) 

we define a site (n x ,n y ) to belong to the A(F)-domain, 
if the Q(n x ,n y ;l) is positive (negative). The parameter / 
grows as a function of the field and so changing the domain 



0.3 




Fig. 9. Magnetoresistance for a random sample of the size N x 
L — 100 x 100 (Ra : Rf : -Rdb = 1 : § : 20). Inset: u(l) defined 
as the density of domain boundaries is here proportional to the 
number of domain boundary resistors in the system for a given 
domain formation. 

distribution, in the very same way as we have seen in our 
very simple model above. Obviously, defining a 9(1) here 
as the number of bonds in the network with a resistor Rdb 
yields a function equivalent to v(l) assumed above (see 
inset of FigEJ. 

The function Q(n x , n y ) is generated from random num- 
bers by smoothening in the same way as the strain distri- 
bution for our mean field discussion. We show the results 
for a simulation for N x L = 100 x 100 and take 



in order to mimic the field dependence of the parameter I 
which should have a qualitatively similar behavior as the 
I in Ea. i|4.3[) . i.e. Z — Z is quadratic in the field for small 
fields and grows faster around the percolation regime (H ex 
dimcnsionless and a = 1/4 and b = 1/2). In Fig|§]we show 
the result for the resistivity with the parameters given in 
the caption. The same qualitative behavior is found here 
as in Fig. [S] including the behavior of v(l), demonstrat- 
ing that the simple phenomcnological model captures the 
basic properties of the role of domain boundaries in the 
transport. 

5 Discussion of the experiment 

Our model reveals two essential points which we can com- 
pare with the experimental situation of Ca,2- x Sr x RuO±. 

(1) We obtained a mean field phase diagram as a func- 
tion of temperature and orthorhombic deformation which 
we can translate to the experimental phase diagram of 
temperature versus Sr-concentration x for 0.2 < x < 0.5. 

(2) We identify our antiferromagnetic phase (with driven 
ferro-orbital correlation) with the region II of the exper- 
imental phase diagram and can discuss some important 
properties. For both points the disorder effect of the alloy 
plays a decisive part. 
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(1) The concentration x = 0.5 is separating region III 
and region II. At this point the system is still dommantly 
tetragonal and the domains with anti-ferro-orbital order 
dominate. The spin susceptibility is Curie-like and shows 
at temperatures ~ IK a irreversible behavior. This low- 
temperature phase was interpreted by Nakatsuji et al. as a 
cluster glass. Within our model we would rather identify 
this phase with an inhomogeneous ferromagnetic phase. 
The observed slow dynamics is then not a glass-like fea- 
ture, but rather an effect due to slow motion of weakly 
pinned domain walls. If the concentration x decreases into 
the region II, gradually the size of the antiferro-orbital do- 
mains is shrinking and the ferromagnetic domains even- 
tually do not percolate anymore, so that the irreversible 
low-temperature phase ceases to exist. 

(2) Close to x = 0.2 in region II the physics is gov- 
erned by the dominance of a percolating domain with 
antiferromagnetic spin and the driven ferro-orbital cor- 
relation. Domains of antiferro-orbital order and the ferro- 
magnetic spin correlation are sparse. In this situation we 
find two important features to compare with experiments. 
The first is the reduced crystal symmetry which due to 
spin-orbit coupling manifests itself in the anisotropy of 
the spin susceptibility. (A puzzle remains in the fact that 
the orthorhombic distortion is untwinncd in all the sam- 
ples investigated so far, so that this anisotropy is very 
pronounced. 0221) The second feature is connected with 
the fact that the size of the minority domains can be in- 
creased by an external magnetic field. As a consequence 
a metamagnetic transition is observed at roughly 2.5T for 
fields along the inplane direction with the maximal spin 
susceptibility and a pronounced peak in the longitudinal 
magnetoresistance around 2.7T. In our previous discus- 
sion we have shown that these features fit well into the 
picture of redistributing domains by the magnetic field. 
In particular, these effects disappear for temperatures ex- 
ceeding 10 K which corresponds to the maximum of the 
uniform spin susceptibility indicating the onset of anti- 
ferromagnetic correlations. Only below this characteristic 
temperature domain formation is possible at all. 

Finally, we would like to comment on the energy scales 
entering our effective Hamiltonian. Our mean field treat- 
ment yields real phase transitions whose onset tempera- 
tures can be taken as the energy/temperature scale. Iden- 
tifying the Need temperature as approximately 10K we 
may choose J to be about 30K. Then, the onset of FM 
order around x — 0.5 is of the order of IK and the critical 
temperature of the AFO order lies between 10 and 20K. 
While the magnetic orders are clearly related to the exper- 
imentally observed behavior of the system, the staggered 
orbital order has not been observed so far. Nevertheless, 
it is important to remark here that the LDA+U calcula- 
tion by Anisimov et al. based on the experimentally deter- 
mined crystal structure data at x = 0.5, clearly reveals the 
staggered orbital correlation as a dominant orbital feature 
including simultaneously ferromagnetism. 16 



6 Conclusions 

The initial motivation to explain the presence of spin 1/2 
observed in the Curie-like susceptibility at the boundary 
between region II and III, has led to the hypothesis of a 
Mott insulator which is restricted to two of three orbitals. 
This would introduce a localized spin 1 /2 and orbital de- 
grees of freedom besides a band of itinerant electrons. 
We have seen that the resulting effective model (Eq.(@J) 
does not only explain the Curie-like susceptibility around 
x = 0.5, but gives a surprisingly good account of a vari- 
ety of additional features: a qualitative description of the 
NM phase diagram of Ca2- ;c Sr. c Ru04 as well as a good 
account of the properties of region II in a magnetic field. 
We do not claim that the features discussed successfully 
by this model is a proof of validity of the initial hypothe- 
sis of orbital selective Mott transition. It is definitely also 
important to take in future studies the coupling between 
the itinerant 7-band and the localized degrees of freedom, 
as it would influence the phase diagram obtained. Never- 
theless, we believe that our discussion provides some evi- 
dence that localized degrees of freedom should be involved 
in the physics of this material. It seems difficult to us to 
explain these characteristic properties of Ca2- :c Sr :r Ru04 
by a model entirely based on itinerant electrons. 
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cussion during the setup of this work. M.T. acknowledges 
the Aspen Center for Physics and the KITP where parts 
of this work have been carried out. We would also like 
to thank for financial support by the Swiss Nationalfonds 
and a grant-in-aid of the Japanese Ministry of Education, 
Science, Culture and Sports. 



A Strain and isospin 

We consider here the symmetry properties of the terms 
containing linear couplings of the isospin to the strain and 
the magnetic field and spin. It is important to notice that 
( )) transform under symmetry operations of the 

Ru02-plane (C^) like the coordinates (x,y). Therefore is 
the basic symmetry property of the isospin components: 

i z = U\+)(+\-\-)H)^x 2 -y 2 

(31) 

7- = |+)(-| + |-)(+| -^xy 

We ignore here the y-component of the isospin, since we 
will not use it. From this follows that the allowed couplings 
to the invariant terms in the coupling to the strain result 
from 61 = e xx — e yy for I z and €2 = e X y for I x . Thus the 
coupling term has the general form 

Kil'ei + K 2 I x e 2 (32) 

with K12 as phenomenological coupling constants. The 
coupling between strain and isospin happens via the crys- 
tal field level splitting. 
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B Spin-orbit coupling 

By the same symmetry argument we can derive the al- 
lowed terms describing the effect of spin-orbit coupling as 
anisotropic coupling of the magnetic field to the spin, 

I z <— * H X S X ~ HySy 

(33) 

J H x Sy HySx 

which leads to 



9l (I z ){H x S x - H y Sy) + g 2 {I x ){H x S y + HyS X ) (34) 

where again gi 2 are phenomenological constants. Note 
that the effect of spin-orbit coupling also involves the itin- 
erant 7-orbital. This is automatically taken care in this 
symmetry consideration. 
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